Stable and unstable regimes in Bose-Fermi mixture with attraction between 

components 
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A collapse of the trapped boson- fermion mixture with the attraction between bosons and fermions 
is investigated in the framework of the effective Hamiltonian for the Bose system. The properties 
of the 87 Rb and 40 K mixture are analyzed quantitatively at T = 0. We find numerically solutions 
of modified Gross- Pitaevskii equation which continuously go from stable to unstable branch. We 
discuss the relation of the onset of collapse with macroscopic properties of the system. A comparison 
with the case of a Bose condensate of atomic 7 Li system is given. 
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I. INTRODUCTION 

Bose-Einstein condensation (BEC) in ultracold atomic 
gas clouds with repulsive and attractive interatomic in- 
teractions 0, 0, [1] have been the subject of intense theo- 
retical and experimental interest in recent years. Besides 
the studies using the bosonic atoms, growing interest is 
focused on the cooling of fermionic atoms [H • Cooling of 
trapped fermionic atoms to a temperature regime where a 
Fermi gas can be considered as degenerate has been pos- 
sible by sympathetic cooling in the presence of a second 
boson or fermion component. Quantum degeneracy was 
first reached with mixtures of bosonic 7 Li and fermionic 
6 Li atoms [H, 0). Later, experiments to cool mixtures 
of 23 Na and 6 Li 0, as well as 87 Rb and 40 K @, i], to 
ultralow temperatures succeeded. The Bose gas, which 
can be cooled evaporatively, is used as a coolant, the 
fermionic system being in thermal equilibrium with the 
cold Bose gas through boson-fermion interaction in the 
region of overlapping of the systems. 

Collisional interaction between bosons and fermions 
greatly affect the properties of the degenerate mixture. 
Theoretical considerations predict the phenomenon of 
component separation for systems with a positive cou- 
pling constant [Toj . instabilities and significant modifi- 
cation of the properties of individual component in a 
case of boson-fermion attraction [ll|, Q3 ■ Effect of den- 
sity fluctuations in a Bose condensate on the fermion- 
fermion interaction with relevant implications for the 
achievement of fermionic superfluidity has been inves- 
tigated in 13[ . The presence of a sufficiently attractive 
boson-fermion interaction can bring about the formation 
of stable fermionic bright solitons fl4j . 

The simultaneous collapse of the two species has been 
observed in a 40 K — 87 Kb mixture by Modugno and co- 
workers Q. K and Rb atoms were prepared in doubly 
polarized states \F = 9/2, rap = 9/2 > and |2,2 >, re- 
spectively. As was shown, as the number of bosons is 
increased there is an instability value Nb c at which a 



discontinues leakage of the bosons and fermions occurs, 
and collapse of boson and fermion clouds is observed. 
The collapse was found for the following critical num- 
bers of bosons and fermions: Nb c ~ 10 5 ; Nk ~ 2 x 10 4 . 
For investigating the stability diagram of the K-Rb sys- 
tem the mean- field model of Gross-Pitaevskii (GP) type, 
coupled to the Thomas-Fermi equation for fermions has 
been used [ill EE EH- The ground state and the sta- 
bility of the system against collapse was considered with 
an imaginary-time evolution scheme. The signature of 
the instability is the failure of the convergence procedure 
toward the ground state of the system, characterizing by 
an indefinite growth of the central density [161 ] . 

In this paper we study the instability and collapses 
of the trapped boson-fermion mixture due to the boson- 
fermion attractive interaction, using the effective Hamil- 
tonian for the Bose system [l7|, [l8[ . The effective Hamil- 
tonian incorporates the three- particle elastic collisions 
induced by the boson- fermion interaction. We ana- 
lyze quantitatively properties of the 87 Rb and 40 K mix- 
ture with an attractive interaction between bosons and 
fermions at T = 0. The stability of this system on the 
basis variational condensate wave function was studied in 
[l8| , and good agreement with experiment [1| was found. 
We estimate the instability boson number Nb c for the 
collapse transition by numerical calculation of the modi- 
fied GP equation and give a comparison with the similar 
picture in a single Bose- condensate with attractive inter- 
action. Our instability analysis will involve dependences 
of the chemical potential /i and the number of boson par- 
ticles N on the value of central density n c of the bose- 
condensate wave function </>(r). Considerations based on 
n c -dependence were introduced earlier in Ref. [l9| for 
studying the stability of a Bose condensate of atomic 7 Li 
in a magnetic trap at nonzero temperature. The calcula- 
tions [l!| confirmed that the gas becomes mechanically 
unstable when the free energy of the system as a function 
of the central density of the gas approaches a maximum 
value. In our case we present arguments based on the 
second variation of the energy functional. We exhibit ex- 



2 



plicitly those variations of the condensate wave function 
S(j}(r) that reveal an instability point of energy functional 
and relate their to a broad class of variations which does 
not preserve the normalization. 



II. EFFECTIVE BOSE HAMILTONIAN 

First of all we briefly discuss the effective boson Hamil- 
tonian [13, EH. Our starting point is the functional- 
integral representation of the grand-canonical partition 
function of the Bose-Fermi mixture. It has the form 
MM: 



Z = J Dl4>*}D[<l>]D[r]Dty}exp{~(SB(cl>*,<t>)+ 

+ s F {r^)+s int {4>*A,r^))} 

and consists of an integration over a complex field <f>{r, r), 
which is periodic on the imaginary-time interval [0, H0\, 
and over the Grassmann field tP{t 1 y) 1 which is antiperi- 
odic on this interval. Therefore, 0(r, r) describes the 
Bose component of the mixture, whereas ip{r,r) corre- 
sponds to the Fermi component. The term describing 
the Bose gas has the form: 



St 



d h 2 V 2 



I'iiW-Blfr.pJ + llfr.r)!'} 



Because the Pauli principle forbids s-wave scattering be- 
tween fermionic atoms in the same hyperfine state, the 
Fermi-gas term can be written in the form: 



S F (r,i>) = J drj dr|^(r,r) (h 
+ y F (r)-/i F )V(r,r)} 



d h 2 V 2 



dr 2m,F 



The term describing the interaction between the two 
components of the Fermi-Bose mixture is: 

S int {r,<t>,r^)=9BF J dr J rfr^(r,r)| 2 |0(r,r)| 2 , 

where gs — AitK 2 a,B /ms and gsF = 2i{K 2 aBF /mi, 
mj = mBrnp I (tub + tuf), tub and mp are the masses 
of bosonic and fermionic atoms respectively, as and olbf 
are the s- wave scattering lengths of boson-boson and 
boson-fermion interactions. K and Rb atoms in the trap 
are in the potentials with an elongated symmetry (A- trap 
asymmetry parameter) 



Vb(t) 



^ {p 2 + Xz% V F (r) = ^( P 2 +Xz>) 

The trap parameters u>b and u> f are chosen in such a way 
that = mpUp/2, so u>f — \/ m B / m f^b ■ Pa- 

rameters hb and /if are the chemical potentials for the 



bose and fermi systems respectively. The chemical poten- 
tial of an ideal fermi gas in a trap is \xp = fkop(6\N f) 1 ^ 3 

The integral over Fermi fields 

Z F = J D[r}D[i>]exp(~(S F (r,il>)+ 

is Gaussian, we can calculate this integral and obtain the 
partition function of the Fermi system as a functional of 
Bose field 4>{t, r) 



Z,- 



exp 



1 



Seff), S eff = / dr / dr/ o// (|0(r,r)|) 



Using the fact that due to the Pauli principle (quantum 
pressure) the radius of the Bose condensate is much less 
than the radius of the Fermi cloud Rp ~ \J Jlf/Vq, one 
can use an expansions in powers of Vp(r)//i F and obtain 
the effective Hamiltonian in the form 



Heff[4>] 



where 



V eff = k a 



dr 



9eff, 



2m B 



|V0| 2 + (K//(r)-MB)|0| 2 + 



leff 



I0| £ 



m B uj% 2 



fc = (l-^Mi/ 2 ffB^)' K 



V2i 



(1) 



3/2 



3ir 2 tf 



9eff 



1/2 2 



9B - ^^Mf 9bf 



2>>c 



9bf 



The first three terms in {T]) have the conventional 
Gross-Pitaevskii [23| form, and the last term is a re- 
sult of boson-fermion interaction. The interaction with 
Fermi gas leads to modification of the trapping potential. 
For the attractive fermion-boson interaction the system 
should behave as if it was confined in a magnetic trapping 
potential with larger frequencies than the actual ones, in 
agreement with experiment [8J. Boson-fermion interac- 
tion also induces the additional attraction between Bose 
atoms which does not depend on the sign of gBF- The 
last term in H e ff (fTJ) corresponds to the three-particle 
elastic collisions induced by the boson-fermion interac- 
tion. In contrast with inelastic 3-body collisions which 
result in the recombination and removing particles from 
the system [24|, |25j, this term for gsF < leads to in- 
crease of the gas density in the center of the trap in order 
to lower the total energy. 



III. NUMERICAL PROCEDURE 

To simplify the formalism we introduce dimensionless 
variables for the spatial coordinate, the energy, and the 



3 



wave function as 

r = a±_r', E = Huj±E', 



where the typical length and energy of the harmonic ex- 
ternal potential are a± = yJh/mBbj^, hu>j_ = Hlob- 

Then the effective Hamiltonian takes the form (the 
primes omitted) 



Heff = J |^|V0| 2 + U 
+ >| 4 + ^| 6 H 



p 2 + Xz 2 



(2) 



where we introduced dimensionless parameters u — 
2g e f f /a 3 L huj± and v = 2g^jFj/a e L huj±. The wave function 
4> is normalized to the number of atoms in the conden- 
sate / d 3 r\cj)(r)\ 2 = N. In the T -> limit considered, N 
coincides with the total number of bosonic atoms in the 
trap. The explicit form of the ground-state wave function 
is obtained by minimizing the energy functional. The 
first order variation of the energy functional gives the 
modified Gross- Pitaevskii equation 



2 



p 2 + z 2 



MB + 7T^r+oH 4 W> = (3) 



The parameters of the 87 Rb and 40 K mixture with 
an attractive interaction between the bosons and the 
fermions are the following |8j: as — 5.25 nm, cibf — 
— 21. 1 nm. The magnetic potential had an elongated 
symmetry, with harmonic oscillation frequencies for Rb 
atoms lo±_ = lub = 2n x 215 Hz and ojb.z — Xojb = 
2ir x 16.3 Hz. At these parameter values characteris- 
tic length a± = 735 nm, chemical potential for fermions 
Pf w ZlfajjB, ujf ~ 1.47 u>b, ko — 1.07, u = 0.11, 
v = —0.0003. Because of a small deviation ko from unity 
we from now on put ko — 1. Note also that we look for 
the ground state of Eq. i.e. the function </>(r) can be 
treated as real one. 

To clarify the main features of instabilities of the sys- 
tem we consider isotropic picture when a problem can be 
considered effectively as a one-dimensional. The case of 
non-spherical symmetry of the trap is recovered at final 
stage by multiplying the critical number of bosons N c by 
the reverse trap asymmetry ratio 1/A. 

So we get the equation (p = ps)- 



A(f) = (r 2 ~2p + u(j> 2 + -W> 4 )^, 



(4) 



Solutions of Eq. (|4|) will be compared to those for the 
single component Bose condensate with attractive inter- 
actions. As an example of Bose system with attractive 
interaction we choose 7 Li [26|. The s-wave scattering 
length is a = — 27.3ao, where ao is the Bohr radius. The 
transverse frequency is lu±/2tt = 163 Hz, so the corre- 
sponding characteristic length is a± = 2.97- 10 -4 cm and 
u = 87ra/a_L = -0.012. 





FIG. 1: Evolution of the profile of the bosonic condensate 
wave function <f>(r) with increasing the central density n' c . 



It is convenient to look for the numeric solutions of Eq. 
introducing the new parameter: central density n c = 
4> 2 (0). Numerical integration of Eq. (j4j with boundary 
conditions 



0'(O) = 0, 



0, 



oo 



defines the family of solutions cf>(r, n c ) depending on cen- 
tral density, the chemical potential p(n c ) being also a 
function of n c . Such an approach when one considers 
n c as an input parameter except for p enables to find 
solutions in the region of instability and to go contin- 
uously from stable to unstable branch in the parame- 
ter space. This approach differs from an imaginary-time 
scheme [27| . where the stability is indicated by requir- 
ing the convergence procedure to the final value. Solving 
Eq. (ID) , one can easily estimate the effective energy E e f f 
corresponding to the functional |2]) and the ground state 
energy E = E e ff + pN, both as functions of central den- 
sity n c . 

Derivative dc/)(r, n c )/dn c is of especial interest, because 
it determines a change in the number of particle through 
the variation in the central density n c 



dN , , 

— = 2 / d A rcp{r,n c ) 
dn r 



d(f)(r, n c ) 
dn r 



If at some n c there comes about dN/ dn c = it results in 
the appearance of zero mode in density fluctuations and 
the onset of instability fl9j | . 

It is convenient to consider <j)(r,n c ) and dcj>(r,n c ) / ' dn c 
as functions of rescaled central density parameter n' c = 
\u\n c . The results are plotted in Fig. Q] and [2] where 
<j>{r, n c ) is supposed to be normalized to 1, and distance r 
is given in units of a± . Figure[T]shows the evolution of the 
profile of the condensate wave function with increasing 
the central density. For comparison, the solution for the 
isotropic harmonic oscillator, 4>h = n~ 3 / 4 exp(— r 2 /2), 
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FIG. 2: Evolution of derivative d<p(r, n c )/dn c at various values 
of central density n' c . Left panel is for the 7 hi system, and 
right panel for BF mixture. 



which corresponds to the ground state of the ideal Bosc 
gas (u = 0, v = 0), is also shown. For n' c < 20, one sees 
the behavior characteristic of the Bose gas with repulsion, 
namely the cloud density becomes more flat at the trap 
center, with increasing radius of the boson cloud. For 
n' c > 20, the solution changes qualitatively: the central 
density begins to increase. Figure [2] shows the evolution 
of the derivative d<f>(r,n c )/dn c for 7 Li-system with at- 
tractive interaction (it < 0) (the left panel) and for BF 
mixture (the right panel). The behavior of BF mixture 
at relatively high densities (n' c > 20) has similar features 
with Lz-system. Then n' c increases there is a continues 
change of the shape of function d<p/dn c - It acquires a neg- 
ative minimum at r < aj_, which results in a saturation 
and a maximum in N(n c ) dependence. 

To relate conditions for stability of a system towards 
small changes in its density profile with thermodynamic 
functions let us consider the total energy E as a func- 
tional of the condensate wave function <p(r) and its gra- 
dient V</>(r) 

E = J d 3 r£(0(r),V0(r)) 

The first order variation 5E should be considered with 
the constraint 



SN = 0, N = d A r\(j){v 



(5) 



As usual we broaden the class of allowable variation using 
the Lagrange procedure with multiplier /i 



E eff = E-fiN, SE, 



eff 







Effective energy density contains one more variable fi 
£ e ff — £e//(<^(r), V<Kr),/i) an d functional E e ff coin- 
cides with the effective Hamiltonian @. At T — func- 
tional E is nothing but the free energy of the system, and 
E e ff is the thermodynamic potential Q = F — fiN. 



Consider functionals E and E e ff taking it on the so- 
lution of Eq. and on particular class of variations 



(r) = <p(r,n c ), 



d<fi(r, n c ) 
dn r 



(6) 



Then they become a function of central density E — 
E(n c ), E e ff = E e ff(n c ). The first order variations 
8E((f>,S(j)) and SE e ff((f>,S(j),ii) considering on functions 
([6]) are nothing but the first derivatives of functions E(n c ) 
and E e ff(n c ). Due to equalities 



SE 



SE, 



eff 



S<f>{r) 5<fi(r) 
there are simple relations [lj 

SE 



= 0, 



d£, 



eff 



dE, 



eff 



dn r 



dE 
dn r 



= d 



= d 



5<p{r) dn c 



dN 

K n e)-r- 

dn c 



SE, 



eff 



S(j){r) dn c 



d£ e ff dfx 
d[i dn c 



d/j, 
dn r 



(7) 



N(n c ) 



which relates extremum points of E(n c ) and E e ff(n c ) 
with an extremum of fJi(n c ) and N(n c ). Note that varia- 
tions (J6j> do not satisfy constraint (f5|), which holds only 
when dN/dn c = 0. It means that © forms a broader 
class of variations and include those of them which do 
not conserve the number of particle. 

Now we relate the behavior of fi(n c ) and N(n c ) with 
the second order variation S 2 E e ff{4>, S<j), VS(f>, fi) taking it 
on functions ([6]). It is related with d 2 E e ff/dn^ through 
the equality 



drA 



S 2 E eff (n c ) + 



+ 



d 2 £ e ff d<f>(r,n c ) dfx dE e ffd 2 \x 



d4>d[i dn c 



dn r 



Taking into account that d 2 £ e f / / 'd(f)dfi 
a simple relation 



S 2 E ef f(n c ) 



dfi(n c ) dN 
dn c dn c 



dfi dn 2 ( 
—2(f> we obtain 

(9) 



Eq. © shows that there is a simple connection between 
S 2 E e ff(n c ) (taken on a particular class of variation) and 
the behavior of ji{n c ) and N(n c ) . At the point of in- 
stability of the system, where dN/dn c = 0, the second 
variation S 2 E e ff{n c ) is equal to zero. As for the second 
variation of functional E, one can write an equality 



S 2 E{n c ) 



dfi(n c ) dN 
dn r dn r 



+ 2fi(n c ) 



^ cU/>(r,n c 



dn r 



d A r 



which involves an additional term. 
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FIG. 3: 7 Li system behavior. Left panels: Chemical potential 
H, number of particle N and ground state energy E as func- 
tions of rescaled central density n' c . Right panels: Chemical 
potential fj,, ground state energy E and effective ground state 
energy E e ff as functions of number of particles N. 



IV. RESULTS AND DISCUSSION 

To compare the qualitative behavior and properties 
near collapse transition of 7 Li- system with BF mixture, 
we have calculated functions fi(n' c ), N(n' c ) and E(n' c ) 
shown at left panels in Fig. 3 and Fig. 4, respectively. 
Right panels show the dependencies fj,, E and E e j f on the 
number of particle. Function fi(n' c ) for 7 Li system (Fig. 
3) gives no sign of singularities near collapse transition. 
The same can be said about E e ff(n' c ) due to Eq. (jSJ. 
However, there is a common feature in behavior 7 Li and 
BF mixture, namely extremums of N(n' c ) and E(n' c ). For 
system with attraction such a point is at n' c0 « 12, and 
for BF mixture at n' cQ « 23. This value of central density 
corresponds to the onset of instability of the system to- 
wards collapse. This feature was recognized in Ref. (l9j 
and connected with the presence of zero mode fluctuation 
of density at this point. 

Note, that the extremum for 7 Li is very wide in n' c . 
Points of extremum in N(n c ) and n(n' c ) is connected 




20 40 60 80 100 2000 4000 6000 



n'c N 

FIG. 4: BF mixture behavior. The left and right panels are 
the same as in Fig. [3] 



due to Eqs. (0 and JSJ with extremum of E(n' c ) and 
E e ff(n' c ). Functions fj,(N), E(N), and E eff (N) have 
forms characteristic of multi-value behavior. Curves with 
asterisks are those parts of /i, E and E e ff which lie 
at n' c < n' c0 and with triangles are those which lie at 
K > n' c0 . 

At n' c < 23 behavior of BF mixture is similar to 
those of Bose gas with repulsion. Numerical results show 
[H, HH that the density profile n(r) can be accurately 
described in the framework of Thomas- Fermi (TF) ap- 
proximation up to n' c ~ 20. In this region the positive 
zero point energy and boson-boson repulsion energy (the 
first two terms in Eq. (JTJ) stabilize the system. How- 
ever, if the central density grows too much, the kinetic 
energy and boson-boson repulsion are no longer able to 
prevent the collapse of the gas. We see similar behav- 
ior of the system with attraction and BF mixture at 
n' c > 23. Likewise the case of Bose condensate with at- 
traction (see, for example, [H, [24], HfJ), the collapse is 
expected to occur when the number of particles in the 
condensate exceeds the critical value Nb c - For BF mix- 
ture curves /i(n), E(N) and E e ff(N) have a point of 
termination which corresponds to the maximum number 
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of particle N cr ~ 6000. Taking into account the asym- 
metry parameter A sa 0.076 we obtain N cr ~ 10 5 which 
is in good agreement with the experiment [26| , A small 
difference in E e f f for stable and unstable branches arises 
solely from a very small difference in chemical potentials 
of these states and not connected with computational 
accuracy. The small difference in chemical potentials of 
this branches (p(N)- curves with asterisks and triangles 
in Fig. 3) is due to a small value of the three- body 
interaction term v. 

To find the disappearance of the local minimum of 
functional E e ff which points to the instability of the 
system, we should explore the second order variation 
S 2 E e ff. 5 2 E e ff changes a sign from positive to negative 
one at the point of instability. In terms of the steepest de- 
scent method absence of local minimum implies that the 
convergence towards the local minimum falls down. The 
second order variation 5 2 E e ff is given by the quadratic 
form on 8<j> and V(5</> and for the functional © has the 
form (for our purposes it is sufficient to consider only real 





FIG. 5: The profile of the condensate wave function </>(r) 
found from the numeric solution of Eq. Q (solid line) and 
the TF approximation Eq. (|10[) (dashed-dotted line). 



eff 



J {{r 2 - 2p + 3u(f> 2 + 5vct) 4 )(S(t>) 2 + (V6<j)) 2 } d 3 r density profile has the form 



Numerical calculations show that S 2 E e ff > on solu- 
tion </> = (j)(r, n c ) if we take 5<fi(r) as a gaussian, satis- 
fying condition ([5]). This implies that an extremum of 
the Hamiltonian is a local minimum. In the case of Bosc 
condensate with attraction the existence of the barrier 
around the metastable state was confirmed in Ref. (30j 
by extensive variational studies of the nearby wave func- 
tion. 

At the stable branch (n' c < n' c0 ) the value of dp/dN is 
negative for 7 Li system and is positive for BF mixture. 
In a homogeneous one-component system (N/V — const) 
dp/dN is proportional to dp/ dp — 1/(p 2 kt) (p = m\(j)\ 2 
is the mass density, >ct is the isothermal compressibil- 
ity of the system) and the criterion of thermodynamic 
stability >ct > reduces to the requirement that dp/ dN 
should be positive. It is easily generalized for an inhomo- 
geneous system which can be treated in the framework of 
local density approximation. In the local density approx- 
imation the density profile n(r) = \4>\ 2 depends on N as 
a parameter and monotonically expands with increasing 
of particle number. So the density n(r, N) undergoes a 
steady increase: dn(r,N)/dN > at any point r within 
a radius of external potential V ex t(r). A density profile 
is determined from the equation 

Pioc(n(r, N)) + V ext {r) = p{N) 

So, in the local density approximation critcrium of 
local stability dpi oc /dn > through the relation 
dp/dN — (dpioc/dn) (dn/dN) gives the stability condi- 
tion dp I dN > for an inhomogeneous system in external 
potential. 

BF mixture at 7 < n' c < 20 safely can be considered 
in the TF approximation (Fig. [29|]. In this case the 



b 2 TF {r) 



T- 



R 2 



R 2 



R 2 = 2p (10) 



where R 2 r = w 2 /(4|?;|), n cr = u/(2|u|), R < R cr and 
is considered at interval < r < R. The evolution of 
the profile corresponds to a monotonic expanding of the 
boson cloud with increasing number of bosons. That is 
why the stable branch of BF mixture corresponds to the 
positive value of dp/dN . 

In contrast, TF approximation is not applicable for 
7 Li system for which the sign of dn(r,N)/dN depends 
on r. At r < a±_ one observes a rapid growth of the 
condensate density, while at r > a±_ a strong deple- 
tion of the condensate density occurs. The negative 
value of dp/dN < in this case can be explained by 
the following. Using Eq. (J9j) the derivative dp/dN in 
terms of the central density can be rewritten in the form 
dp/dN = (dp/dn c ) {dn c /dN) = S 2 E eff {n c )/(dN/dn c ) 2 . 
So the value of dp/dN gives, at most, an information 
about a particular class of variations © of the functional 
@, which do not satisfy ([5]). For variations satisfying 
condition ([5]) we have 5 2 E e ff > 0. So solution (f>(r,n c ) 
provides the local minimum of E e ff on those functions 
which preserve the normalization. Stability of the system 
is explained by dynamical reasons, namely for collapse of 
system to occur, fluctuations with SN ^ should exist, 
which energetically are not favorable at temperatures un- 
der discussion. 

In conclusion, our analysis of the stability of K-RB 
Fermi-Bose mixture on the basis of effective Bose Hamil- 
tonian shows the clear resemblance to the behavior of 
7 Li system. There is a value of central density at which 
small variations of density profile conserve the number of 
particle SN = and the second variation 5 2 E e ff changes 
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the sign. The value we determine for N c is in very good 
accordance with experiment. Points of extremum of func- 
tions /z(n c ) and N(n c ) is related with the first and the 
second derivatives of functions E(n c ) and E e ff(n c ). 

We note that the investigation of the actual dynamics 
after the system has been driven into the unstable region 
would require a description that go beyond the stationary 
scenario of Eq. (|U), in similar fashion to what happens 
during the collapse of a single Bose-Einstein condensate 
with attractive interaction [24], [25| . Here we will not dis- 
cussed these aspects, since we are concerned with the 
determination of the critical values for the onset of insta- 



bility. Another interesting issue concerns the relevance 
of finite temperature effects, which are not included in 
the present treatment. 
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